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Abstract 

A growing or compressed thin elastic sheet adhered to a rigid substrate can 
exhibit a buckling instability, forming an inward hump. Our study shows that 
the strip morphology depends on the delicate balance between the compression 
energy and the bending energy. We hnd that this instability is a first order 
phase transition between the adhered solution and the buckled solution whose 
main control parameter is related to the sheet stretchability. In the nearly- 
unstretchable regime we provide an analytic expression for the critical threshold. 
Compressibility is the key assumption which allows us to resolve the apparent 
paradox of an unbounded pressure exerted on the external wall by a confined 
flexible loop. 


1 Introduction 

Elastic strips and rods are fundamental physical structures that can be found in 
many contexts on many different length-scales. Often, they appear in the the¬ 
oretical description of physical phenomena where the interplay between physics 
and geometry plays a key role. The collapse of an elastic sheet adhered to a 
curved substrate is certainly one of them. This discontinuous buckling insta¬ 
bility appears in various problems in biology (shapes and morphogenesis of the 
mitochondrion (Sakashita et ah, 2014)), physics (packing of thin hlms (Bone 
et ah, 2006), delamination of thin elastic films (Vella et ah, 2009; Wagner and 
Vella, 2013), design of stretchable electronic devices (Sun et ah, 2006), elasto- 
capillary snapping (Fargette et ah, 2014)) and engineering (collapse of buried 
steel pipelines (Omara et ah, 1997; Vasilikis and Karamanos, 2009)). 

A simple example illustrates the physical phenomenon we are concerned 
with. Let us imagine to insert a rectangular piece of paper inside a rigid cylin¬ 
drical substrate. We hold the two edges with the hands so that the sheet of 
paper is kept in contact with the cylinder. Subsequently, we apply an increasing 
tangential compressive force. Due to presence of the surrounding container, the 
sheet of paper cannot freely deform to the outside. Initially it appears to be 
undeformed or slightly compressed. However, when the applied force crosses a 
suitable threshold the strip abruptly buckles and forms an inward hump. What 
is the critical force at which the collapse occurs? How are the geometric features 
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of the hump at the transition related to the material constants and to the ge¬ 
ometry of the container? Although the buckling of slender elastic structures is 
one of the oldest fundamental problems, nobody seems to have addressed these 
points. Furthermore, it is worth noticing that classical bifurcation analysis fails 
in this problem due to the presence of the conhnement. 

A strictly related problem deals with the packing of a flexible cylinder in¬ 
side a rigid circular tube of smaller radius. This is widely explored in the 
literature, both theoretically (Cerda and Mahadevan, 2005) and experimentally 
(Boue et ah, 2006). Recent papers have also included the possibility of adher¬ 
ence by capillarity (De Pascalis et al., 2014; Rim et ah, 2014) or the extension to 
the spherical geometry (Kahraman et ah, 2012). In the usual theoretical treat¬ 
ment, the strip is modelled as an inextensible Euler beam, whose energy is only 
associated to the bending mode. The inextensibility assumption, first suggested 
by Rayleigh (1877), is based on energetic considerations. The bending energy 
and the stretching energy scale differently with the beam thickness. When the 
structure thickness is diminished without limit, stretching is energetically pro¬ 
hibitively expensive and the distortion is pure bending. However, the validity of 
this assumption is challenged when the thickness becomes comparable to some 
other length- scales in the problem. Moreover, axial compressibility is known to 
change the nature of the phase transition. For example, the usual supercritical 
bifurcation point in the classical inextensible Euler beam becomes a subcritical 
bifurcation (Magnusson et ah, 2001; Audoly and Pomeau, 2010). 

However in our confined problem, the assumption of an axially inextensible 
beam shows its inadequacy in a more fundamental way. For instance, in the 
limit of arbitrarily small hump, the inextensible Euler’s Elastica model predicts 
the unphysical result of an inhnite pressure exerted by the beam on the external 
container (Cerda and Mahadevan, 2005; De Pascalis et ah, 2014). This prevents 
the transition to a buckled state as it would require an inhnite force to induce 
the buckling. 

To avoid this unacceptable behaviour, in this paper we revise this effect 
by allowing the sheet to both stretch and bend. Within this framework, the 
energy landscape exhibits a hrst order phase transition between a state where 
compression energy dominates, to another where bending energy prevails. From 
the mechanical point of view, a snap through buckling occurs: as the sheet grows, 
it suddenly switches from a completely adhered solution to a buckled solution 
showing a symmetric inward hump. The transition is studied numerically and, 
in the limit of weak stretchability, we provide the analytic expressions for the 
main quantities at the critical point. Interestingly, we find that for homogeneous 
beams the threshold does not depend on the material parameters but only on 
the geometric features of the problem. 


2 The model 

We assume a planar configuration so that the longitudinal prohle of the strip 
can be modelled as a stretchable and flexible rod. This is represented by a 
parametric curve r(S'), with S € [—1//2,1//2], where L denotes the length of the 
strip in the stress-free configuration and S' is the referential arclength. The gov¬ 
erning nonlinear Kirchhoff equations express the balance of linear and angular 
momentum. When an external distributed force f is applied and in absence of 
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external distributed torques, they read 


T' + f = 0, (1) 

M' + r' X T = 0, (2) 

where T(5') and M(S') represent the internal forces and the internal torques, 
respectively. The prime denotes differentiation with respect to S; hence, r' = 
At, where A is the local stretch and t is the unit tangent. The inextensibility 
condition corresponds to A = 1. 

Balance equations are complemented by the linear constitutive laws (Dill, 
1992): 


T — 6(A — l)t + Tnii, 

(3) 

IS/I = kt X t', 

(4) 


where the positive constants b and k represent the stretching and the bending 
rigidity, respectively. The quantity £ = \/k/b defines an intrinsic characteristic 
length. Since b = EA and k = El, where E is the elastic modulus of the 
material, A is the area of the film cross-section and I is the second moment 
of area of the strip cross-section, the intrinsic length I turns out to be of the 
order of the strip thickness. Only the tension is given constitutively; the shear 
internal force T„ is related to the derivative of the internal torque by Eq. (2). 
We parametrize the tangent and the normal unit vectors by 

t = cos6{S)ex + sm8{S)ey, n = — sin 6 '(S)e 2 , -|- cosd{S)ey 

and, hence, = t x n. Consequently, Eqn. (4) reduces to M = kO'e^. 

Furthermore, we posit a symmetric equilibrium configuration of the rod and 
assume that the displacement of its end-points can be controlled: the point with 
referential coordinate S = L/2 undergoes a tangential compressive displacement 
rAa, Atr > 0. The strip can both stretch and bend, therefore there are at 
least two equilibrium configurations with the same imposed Aa: (i) the adhered 
solution where the strip, uniformly compressed, totally adheres to the wall, and 
(ii) the collapsed solution with an inward symmetric hump. The adhered solution 
is a uniformly compressed arc of circumference with 

^adh 1 2rAcyjL, ^adh — ^adh jr~ 

On the other hand, the collapsed solution comprises two parts: a free (non- 
adhered) curve for S € {—S,S) and two adhered pieces for S € [—L/2,—S] U 
[S, L/2] (see Figure 1). The symmetry of the strip implies that the function 
0{S) is odd and allows us to restrict the study of the solution in the range 
Se[0,L/2]. 

Since the adhered part, S € [S', L/2], is an arc of circumference with 8'(S) = 
—A/r, the internal torque is M = —k{\/r)ez. Therefore, the equilibrium equa¬ 
tions for this imposed geometry necessarily give T„ = 0 and = 0, whence 
A(S) = A = const. However, the compressive tangential force F = —Ft for a 
given Aa is not known a-priori but it is related to the stretch by F = fe(l — A). 
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Figure 1: Schematic representation of the strip deformation. The gray curve repre¬ 
sents the referential configuration, while the blue curve shows the buckled configura¬ 
tion with a hump due to the displacement rAa of the end points. 


In the non-adhered region, S € [0,5], there is no applied distributed load 
(f = 0). Furthermore, since 0(S} is odd, the vertical component of the tension 
vanishes, while the horizontal component, denoted by h, is constant and satisfies 
the identity b{\ — 1) = h cos 6. The balance of the torque yields the second order 
differential equation for the free part 

9” + Xt sin 6 = 0. (5) 


where r = —h/k. 

The deflection angle 6{S) should satisfy the Dirichlet boundary conditions 
0(0) = 0, 6{S) = 6. Further boundary conditions are provided by requiring the 
continuity of the tangential component of the internal force and the continuity 
of the torque at the detachment point S = S. These yield 

/icos0 = —6(1 — A), 9'{S) = —Xr~^, 

respectively. One more constraint translates the geometrical condition that the 
detachment point must lie on the circumference of radius r 

A cos 0d5 = —r sin 0. (6) 

This establishes a relationships between S and 0. The displacement of the end¬ 
point, rAa, is related to A, S and 0 by (Figure 1) 

rAa = r0-h(l-A)^-tA5. (7) 

This geometrical identity is simply obtained by requiring that the total length 
is the sum of the hump length and the adhered length. 
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3 Results 


We start by observing that for a perfectly unstretchable strip the compression 
modulus diverges and the intrinsic length vanishes, i! = 0. The tangential com¬ 
ponent of the internal force, Tt, then becomes a Lagrange multiplier associated 
to the inextensibility constraint A = 1 and the adhered solution is admissible 
only for Aa = 0. 

By contrast, when both compression and bending are allowed, the energy 
must attain its absolute minimum in the observed solution, either adhered or 
collapsed, which is then expected to be stable. The other possible solution is 
either non-existing, unstable or metastable, depending on the energy landscape. 
The total elastic energy consists of two terms, related respectively to the bending 
mode and to the stretching (or compression) mode: 

k /■^/2 b 

w=- {e'fds+- ii-\fds. ( 8 ) 

^ J-L/2 ^ J-L/2 

Equation (8) assumes a particularly simple analytic form when evaluated in the 
adhered solution 

fE.., = ^^ + 2^Aa[-l + ^(l + r>«], (9) 

where ^ = i/r > 0 measures the compressibility strength at fixed radius. 



Figure 2: Plot of the equilibrium free energy W for the adhered solution (red) and the 
collapsed solution (blue) with L — |7rr and ^ = 0.01. Solid lines represent solutions 
with least energy; dashed lines represent metastable states or maxima. For small Aa 
only the adhered solution is admissible. The dotted line represents the energy of a 
perfectly inextensible strip as given in Eq.(17). 

On the other hand, the energy associated with the collapsed solution, Wcoi, is 
not amenable for a simple analytical approximation and is best calculated with 
a numerical simulation. Figure 2 shows the comparison between the energies 
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associated with each branch as a function of Aa, for fixed ^ = 0.01. The adhered 
branch depends quadratically on Aq (see equation (9)). The collapsed solution 
comprises two branches. The upper branch, with higher energy, corresponds to 
small humps, while the lower branch corresponds to equilibrium solutions with 
larger humps. Finally, it is worth noticing that for sufficiently small values of 
Aa only the adhered solution is admissible. 

Figure 2 clearly shows the first-order transition that occurs between the two 
solutions. For small displacements of the end-points, the adhered solution is 
energetically favoured. However, when the displacement crosses a critical value 
r{Aa)cr, the energy attains its absolute minimum at the collapsed solution with 
larger hump. Therefore, as Aa is increased, the strip undergoes a discontinuous 
transition, abruptly passing from an adhered configuration to a buckled solution 
with an inward hump. The presence of compressibility implies a lower bound for 
the hump size: humps whose dimensions are small compared to the characteristic 
length £ cannot be observed and the strip is simply compressed. The critical 
threshold increases as the strip becomes softer (see Figure 3). In particular, in 
the inextensible limit the critical threshold vanishes, as expected. 



Figure 3: Critical threshold (Aa)cr as a function of the reduced compressibility. 
The solid line represents the numerical results, while the dotted line represents the 
analytical approximate result as given in Eq.(18). Typical shapes are shown in the 
inset. 


In the nearly unstretchable regime (^ ^ 1) one can estimate the critical 
threshold analytically. Indeed, at low ^ and close to the critical point, it has 
been observed numerically that the collapsed solution is nearly equivalent to 
the buckled profile of a perfectly inextensible strip. Figure 2 clearly shows this 
agreement in terms of elastic energies. From the analytic approximation we also 
learn that the relevant variable for scaling the critical threshold is the reduced 
compressibility, defined as 

Co := 

This indicates an unexpected non-trivial dependence on the geometrical and 
material parameters. For homogeneous materials is just a geometrical di- 
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mensionless parameter, since both the compression and the bending rigidity are 
linear functions of the Young modulus. By contrast, the material moduli could 
play a role in the case of composite materials where the dependence of b and k 
on the Young modulus is more complex. 

When the compression modulus Is large <C 1), compression becomes ener¬ 
getically prohibitively expensive for small Aa. The strip then prefers to buckle 
and pay some bending energy provided it can relax its compression energy. 
Therefore, it is natural to approximate the collapsed strip as inextensible. The 
advantage of this approximation is that it is now possible to perform an ana¬ 
lytic treatment. In particular, Eq.(5) simply reduces to the nonlinear pendulum 
equation with first integral 

{d')^ = 2 t{cos9 — cos6o), (10) 


where 6o = 9{So) € [0, rr] is the maximum value of 9{S) in (0, S). The boundary 
condition for 6' at S becomes 9'{S) = —Ijr. This is used together with the 
first integral to eliminate r in favour of Oq. Thus, we arrive at the first order 
differential equation 


9 ' = ±- 


1 / cos 9 — cos 9o 


r V cos 9 — cos 9o 


( 11 ) 


where the sign + (respectively, —) is to be used in the interval S € (0, S'o) 
(respectively, S € (5'o,S)). By symmetry 0(0) = 0 and Eq.(ll) evaluated at 
5 = 0 shows that cos 9 — cos0o > 0. This gives a restriction on the possible 
values of 9: |0| < 9o. Furthermore, Eq.(ll) is an ordinary differential equation 
which can be solved by separation of variables in (0, S) (see (De Pascalis et ah, 
2014, Sect. 3)) 


2F(go) - F{q) 


S I 1 — cos 00 
2r V cos 0 — cos 00 ’ 


( 12 ) 


where F denotes the incomplete elliptic integral of first kind (Abramowitz and 
Stegun, 1970) and, for ease of notation, we set 


qo ■- {6'o/2,csc^(0o/2)}, q := {0/2,csc^(0o/2)}. 


Similarly, we simplify Eq.(6) (with A = 1) as follows 


2E(go) - E(g) = 


S cos 9o + r sin 0 j 1 — cos 0o 
2r (1 — cos 0o) V cos 0 — cos 0o ’ 


(13) 


where E represents the incomplete elliptic integral of second kind. Equations 
(12), (13), and the inextensible version of Eq.(7), rAa = r9 + S, provide S, 
00 and 0 as functions of Aa. In the limit Aa <C 1, after some tedious but 
straightforward calculations which we omit for brevity, it is possible to find the 
asymptotic expansions of 0o and 0 (De Pascalis et ah, 2014) 


00 « 2.3454(Aa)^/® - 0.78762Aa -b 0.68361(Aa)®/®, (14) 

0 « -2.2894(Aa)^/® -b 0.83072Aa - 0.62489(Aa)®/^ (15) 


The energy of an inextensible buckled solution has a particularly simple expres¬ 
sion when it is written in terms of 0o, 0 ans S. To this end, we insert Eq.(ll), 
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Eq.(6) and A = 1 into Eq.(8) to get 


Wcoi = k 


{e'fdS + k 


r-i/2 


Jo 

2 


(sVds.i 


cosO{S) — cos 00 
cos 6 — cos 00 


k S cos 0 + r sin 0 
cos 0 — cos 00 


dS+X(^-S 


(16) 


The asymptotic expansions of 0o and 0 are then inserted into Eq.(16) to yield 
the approximate expression for the energy of the collapsed solution 

Wcoi~^^ + ^{wo{Aay^^ + wiAa + W2(Aa)®^®), (17) 

where wq « 23.113, Wi « —18.532, W 2 ~ 1.6131. Thus, we are able to give an 
approximate analytic solution of the transition equation Wcoi = Wadh, which to 
leading order reads 

(Aa)cr. « 6.5813 (18) 

As shown in Figure 3, this expression provides a very good account of the nu¬ 
merical results, in the range of compressibility considered. 

Figure 4 sketches the critical behaviour of the detachment angle. This angle 
is directly related to the hump dimension and is easily accessible experimentally. 
For a given |0cr| represents the minimum stable detachment angle. Lower val¬ 
ues of 0 corresponds to metastable or unstable solution. Thus, we can conclude 
that there is a limiting size below which the hump cannot develop. The dotted 
line in Figure 4 represents the asymptotic approximation \6cr\ ~ 4.29 ■ 

The lower bound for the blister size corresponds to a higher bound for the 
applied tangential force and thus for the pressure exerted on the delimiting 
container. Figure 5 reports the behaviour at the transition of the tangential 
force for the adhered (red line) and collapsed (blue line) solutions, as functions 
of the reduced compressibility. Since at the transition the solution passes from 
one to the other, the tangential force undergoes a discontinuous jump. This 
is consistent with the idea of buckling as the mechanism through which the 
system relaxes its internal stress. For fixed ^o, the correspondent value of the 
red line shows the maximum internal force acceptable by the system under 
pure compression. The critical load at which the buckling occurs is given by 
E’er = 2br{Aa)cr/L- Furthermore, when Aa is below its critical threshold the 
internal force is simply proportional to the strain and hence it vanishes as Aa 
tends to zero. 


4 Concluding remarks 

In many applied contexts it is useful to study the slightly different problem 
where the sheet Is closed, i.e. the lateral ends of the strip are glued together. 
For instance, the elastic sheet can be made of a growing biological membrane 
confined by a rigid cylinder of radius r. In this case, the circumferential growth, 
suitably measured by e ~ {L — 27rr)/(27rr), may trigger the instability. For 
moderate growth, the excess material can simply result in a uniform compression 
of the adhered strip and in an increase of the hoop stress. By contrast, a further 
increase of the growth (when the critical threshold has been crossed) leads to 
the buckling of the membrane and the sheet is then only in partial contact 
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Figure 4: Detachment angle at the transition as a function of The dotted line 
represents the asymptotic approximation. 


with the container. Our results can be simply adapted to this closed problem. 
As expected we again find that there is a critical growth Ccr above which the 
delamination is preferred with respect to the adhered solution. Indeed, the 
equations governing the post-buckling behaviour are the same except for Eq.(7) 
that is to be replaced with 

■Krt = re+{l-X)L/2 + \S. 

It is worth noticing explicitly that even the energy expression (8) remains un¬ 
changed. By contrast, the adhered solution is characterized by a uniform stretch 
Xadh = 1/(1 -f e) and its elastic energy is 

Wadh = 7rfcr“^(l -I- e)“^(l -|- 

However, it turns out that the two problems are mathematically equivalent 
provided we substitute Aa with 7re. We finally find the following approximation 
for the critical growth: ter ~ 41.0912 

As an application of our results, we look at the experiment on the packing 
of flexible films reported in Boue et al. (2006). In the first part of their paper, 
the Authors studied the pressure exerted by the sheet on the external container 
as a function of e, here named the confinement parameter. Their experimental 
data exhibit an increasing pressure as the confinement parameter t is decreased. 
Numerical studies of the elastica model indicate indeed a divergent pressure in 
the limit of e —> 0. This is also consistent with the theoretical results (Cerda 
and Mahadevan, 2005; De Pascalis et ah, 2014), which in addition show an 
asymptotic dependence of the internal forces, and hence the pressure, on t~^^^. 
However, the thickness of the sheet used in the experiment is h — 0.1mm, while 
the container radius is r = 26 mm. Since for a rectangular section £ = hf\/l2, 
we obtain that ^ = 1.11 x 10“^ and ter = 1-17 x 10“^. Unfortunately, the 
smallest experimental value of e reported in Figure 3 of Boue et al. (2006) is 
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Figure 5: Tangential compressive force at the transition in the adhered solution (red 
line) and collapsed solution (blue line) as a function of The tangential force has a 
finite jump at the transition. The dotted lines show the asymptotic approximations. 


just above this critical value and hence they are not able to see the transition 
experimentally. On the one hand, this justify the validity of the Elastica model 
above the critical threshold. On the other hand, our results predict that, by 
further decreasing e, the hump should suddenly disappear. At the transition, the 
pressure (per unit length) on the container undergoes a jump to Per ~ br~^eer, 
while below €cr the pressure is proportional to e, hence it goes to zero as e —>■ 0. 
This solve the apparent paradox, described in Cerda and Mahadevan (2005) , 
De Pascalis et al. (2014) and Boue et al. (2006), of a diverging pressure in the 
low confinement regime. 

In conclusion, the present analysis recognizes the key role played by stretch- 
ability for a correct description of the mechanical collapse of a loaded or growing 
elastic thin sheet adhered to a curved substrate. Even when the sheet is nearly 
unstretchable, the energetic contribution of compression cannot be neglected 
when the length scales in the problem are comparable with the sheet thickness 
(i.e., small humps). Our model yields analytic approximations for the critical 
threshold and the minimum size for the humps. This threshold depends on the 
reduced compressibility that, for homogeneous materials, is a purely geometrical 
dimensionless parameter. Furthermore, it sets an upper bound to the exerted 
pressure on the container. 

Future works can consider the effects of other physical quantities (i.e. cap¬ 
illary adhesion, intrinsic curvature and gravity) or the coupling of growth with 
internal stress (Tiero and Tomassetti, 2014) on the morphology and the critical 
parameters. 
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